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Abstract A quasi-simple shear test, which is the most direct 
method for examining the shear properties of sheet metals, 
has been applied to measure the shear moduli of wood. 
Buna (Feigns crenata Blume) with variously sized shear re¬ 
gions was used for the test specimens. Strain gauges were 
mounted in the center of the shear regions to measure the 
shear strains. The shear tests were carried out to determine 
the shear moduli in the radial and tangential planes. Appar¬ 
ent shear moduli obtained from the experimental results 
were corrected by finite clement method (FEM) simulation 
of the shear region, where both shearing and bending are 
produced. It was found that the corrected shear moduli are 
roughly independent of test conditions, and their values 
are in good agreement with the data obtained from 
bending-shear tests. This suggests that the method em¬ 
ployed here can effectively estimate the shear moduli of 
wood. 

Key words Shear modulus ■ Quasi-simple shear test ■ Ap¬ 
parent shear modulus • Stress analysis 


Introduction 

The torsion test and off-axis tension or compression test are 
commonly used to measure the shear modulus of wood, 1 
although in some cases the static bending test or flexural 
vibration test is utilized instead. 2 ' 3 Research on application 
of the Iosipescu shear test to this measurement was also 
reported. 4 Having different features, these conventional 
tests are used depending on the requirements for accuracy 
and simplicity of the experiment. 
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In this paper, the feasibility of a simple shear test as an 
alternative method for evaluating the shear modulus of 
wood is discussed from the aspects of both experimental 
and finite element method (FEM) simulation. Although the 
simple shear test is basically the most direct method for 
investigating the shear properties of materials, it is difficult 
to obtain the shear modulus of wood by this method. This is 
because we cannot measure exactly the small deformations 
in the shear regions of the test specimen without accurate 
equipment. To overcome the experimental difficulty, strain 
gauges were used here. In general, shearing and bending, 
which is a minor deformation, take place simultaneously in 
the shear regions. The more the shear regions are governed 
by shear, the more fully are the simple shear deformations 
attained. The possibility of simple shear depends on the 
dimensions of the shear zones. Taking this point into ac¬ 
count, we propose a simplified procedure to obtain the 
proper shear modulus value. 


Experiment 

Specimens 

Buna (Feigns crenata Blume) was used for the specimens. 
The density and moisture content of the specimens were in 
the range of 0.71-0.76g/cm 3 and 6.9%-8.7%, respectively. 

Quasi-simple shear tests 

The simple shear test was originally developed for sheet 
metals. 5 The features of this test as shown in Fig. 1 are that 
shear zones are disposed symmetrically to produce a stable 
deformation and that only simple shear is expected to occur 
in the shear zones with the strong restraint as the specimen 
dimension c becomes sufficiently small. 

The test specimen used here is shown in Fig. 2. As illus¬ 
trated in Fig. 3, the specimen, having a thickness of 5 mm, 
was fixed tightly to the shearing test jigs with eight bolts to 
generate the shear deformations within the shear zones, 
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Fig. 2. Geometry of test specimen with thickness t of 5 mm (c = 10, 
15 mm) 

where c is 10 and 15 mm. In this case, because the coexisting 
bending cannot be ignored, we defined this test as the quasi¬ 
simple shear test. To measure the shear strains in the center 
of the shear zones, biaxial strain gauges (gauge length 2 mm, 
base circle diameter 7mm) (Tokyo Sokki) were bonded 
on the specimen’s surface. The corners of the shear zones 
were rounded to a radius of 1mm to avoid the stress 
concentrations. 

We measured the shear moduli in the longitudinal-radial 
(LR) and longitudinal-tangential (LT) planes. Two speci¬ 
mens were prepared for each measurement: One test speci¬ 
men was loaded in a direction parallel to the longitudinal 
direction (LR specimen), and a second one was loaded 


perpendicular to the longitudinal direction (RL specimen). 
Great care was taken during specimen preparation to 
achieve a radial or tangential plane as precisely as possible, 
particularly in the LT specimens. The tests were performed 
by applying the shearing load in a computer-controlled test¬ 
ing machine (Shimadzu AG-10TA) at a constant crosshead 
speed of 0.05mm/min. The measured values of load and 
strains were then sent to a personal computer via a data 
logger (NEC San-ei D5100). Five specimens were used for 
each test condition. 

Bending-shear tests 

Four-point bending tests were also carried out for compari¬ 
son purposes. The testing method is shown in Fig. 4. Beam 
specimens with the dimensions of 25 mm (radial direction) 
X 25 mm (tangential direction) X 480mm (longitudinal 
direction) were prepared from the same lumber. To mea¬ 
sure the shear strains in the middle between the loading 
point and support, biaxial strain gauges, which are the same 
as those used in the shear tests, were bonded on the speci¬ 
men. Loading heads with a radius of 30 mm were used. The 
speed of loading was 1 mm/min. Five specimens were used 
for each test condition. 

FEM simulation of quasi-simple shear 

Quasi-simple shear tests were simulated by a three- 
dimensional FEM program 6 modified for orthotropic mate¬ 
rials. Figure 5 shows the FEM model for the LR or RL 
specimens with c at 10 mm, where the z-axis is taken in the 
loading direction. The models were divided into the finite 
elements having the dimensions of 5 X 1 X 1mm. The 
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Strain gauge 



Fig. 5. Coordinate system and finite element mesh for shear zone of 
LR/RL specimen. L, R, T, longitudinal, radial, and tangential direc¬ 
tions, respectively 

lowermost elements lying outside the shear zone in the y- 
axis were fixed, whereas the uppermost ones were given a 
prescribed displacement X along the z-axis. The elastic con¬ 
stants used for RL specimens are shown in Table 1. These 
values were obtained from the literature 7 except for 
Young’s modulus E y , which was obtained from compression 
tests. The listed shear modulus G yz , which is the most effec¬ 
tive parameter on the shear deformations, is regarded as the 
average value for buna. 


Results and discussion 

Numerical stress distribution 

Figure 6 shows the distributions of the shear stress r yz for X 
of 0.1mm in the LR plane. Compared with the LR speci¬ 


fic 6. Distributions of shear stress. Stress values were obtained at 
points x = 0, y = c/2 + 1 mm when the prescribed shear displacement 
l is 0.1 mm 


Table 1. Elastic constants used for FEM calculations 


Parameter 

Constant 

Young’s modulus (GPa) 

E x 

0.588 

E, 

11.760 

E- 

1.176 

Shear modulus (GPa) 

G„ 

0.637 

c y! 

0.980 

G« 

0.196 

Poisson’s ratio 

V 

y xy 

0.025 

V 

yz 

0.400 

V 

y ZX 

0.650 


The data are for a specimen loaded in a direction perpendicular to the 
longitudinal direction (RL) 

FEM, finite element method 


mens, the RL specimens have the approximately constant 
values of x yz over a wide range of the z-axis, and these values 
are nearly equal to (r vz ) simple , which denotes the shear stress 
for simple shear and is given by 

(rJ = G yz ■ — (1) 

V /simple y c 

For the LR specimens, the shear stresses in the central 
part of the shear zone on which the strain gauge is mounted 
are also close to (r y .) simple ; the other stress components can 
be neglected there. Similar tendencies were observed for 
the shear stresses in the LT plane. Figure 7 shows the rela¬ 
tion between the ratio of the average shear stress (r vz ) av to 
( r vz)simpie and the dimension c. The stress (r yz ) m was defined 
as 

Mav = (2) 

where F, is the summation of the nodal forces associated 
with the specified displacement X. As mentioned above, the 
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c (mm) 

Fig. 7. Influence of shear zone size c on the degree of simple shear 
deformation 

( r y;)av/( T 'y Z )simpie ratio, which is unity for simple shear, de¬ 
creases linearly with increasing c because of bending. 

Apparent stress-strain relation 

Figure 8 shows an example of apparent shear stress-strain 
relations obtained from the quasi-simple shear tests. The 
apparent shear stress (r y .) exp and the shear strain y yz were 
given by 


( T ^)exp 

P 

2Tt 

(3) 


>2 - ft) + («3 - £ 4 )] 

(4) 


where P is the applied load, and £,-£4 are the normal strains 
measured here. The apparent shear modulus (G yz ) sxp was 
determined from the initial linear segment of the apparent 
shear stress-shear strain diagram. 

Correction of apparent shear modulus 

The shear modulus must be evaluated with the actual shear 
stress (r v .) ctr in the central part of shear zones. Here we 
assumed that the apparent shear stress (r vz ) exp is equivalent 
to the average shear stress (r vz ) av obtained from the numeri¬ 
cal results by FEM, and that (r yz ) clr is approximated by the 
calculated value. According to these assumptions, the ap¬ 
parent shear modulus ( G yz ) exp was corrected to estimate the 
proper value of the shear modulus (G yz ) est as follows: 


( G -)est “ a '( G «)exp 

(5) 

a = ( r *L/M„ 

(6) 


Because (G yz ) est obtained from Eqs. (5) and ( 6 ) generally 
differs from the provisional value (Table 1) used for FEM 
calculations, an iterative approach was employed. Once the 
shear modulus (G yz ) est is obtained, it is used to obtain a more 
proper value of a in the next step (Fig. 9). It was found that 



Fig. 8. Typical relations between apparent shear stress and shear strain 



Fig. 9. Flow chart of determination technique of correction factor a c 
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Table 2. Change of correction factor a c with Young’s modulus (c = 10 mm) 


Specimen 

Young’s 

modulus (GPa) 

Shear modulus (GPa) 

a c 


E x 

E, 

E z 

G xy 


G„ 

LR 

0.588 

0.882 

11.76 

0.196 

0.98 

0.637 

1.52 



1.176 





1.45 



1.47 





1.41 

RL 

0.588 

8.82 

1.176 

0.637 

0.98 

0.196 

1.12 



11.76 





1.10 



14.7 





1.09 

LT 

1.176 

0.392 

11.76 

0.196 

0.637 

0.98 

1.55 



0.588 





1.50 



0.784 





1.46 

TL 

1.176 

8.82 

0.588 

0.98 

0.637 

0.196 

1.10 



11.76 





1.08 



14.7 





1.07 


LR, loaded in a direction parallel to the longitudinal direction in the radial plane; LT, loaded in 
a direction parallel to the longitudinal direction in the tangential plane; TL, loaded in a direction 
perpendicular to the longitudinal direction in the tangential plane 


Table 3. Apparent and corrected shear moduli for various test 
conditions 

Specimen 

c 

(mm) 

(^yz)exp 

(GPa) 

CV 

(%) 

a c 

(^yz)est 

(GPa) 

LR 

10 

0.89 

6.7 

1.45 

1.29 


15 

0.75 

5.4 

1.54 

1.16 

RL 

10 

1.10 

3.6 

1.10 

1.21 


15 

1.04 

4.7 

1.18 

1.23 

LT 

10 

0.63 

3.5 

1.50 

0.95 


15 

0.61 

3.9 

1.55 

0.94 

TL 

10 

0.84 

5.3 

1.08 

0.91 


15 

0.81 

4.5 

1.14 

0.93 


Table 4. Apparent and corrected shear moduli obtained from the four- 
point bending test 


Specimen 

S~< 

'- 7 exp 

(GPa) 

CV 

(%) 

a N 

a F 

s-y N 
'- 7 est 

(GPa) 

Gl, 

(GPa) 

RL 

0.82 

6.3 

1.60 

1.50 

1.31 

1.22 

TL 

0.59 

3.9 

1.56 

1.50 

0.92 

0.89 


Superscripts N and F indicate the author’s method utilizing FEM simu¬ 
lation and Fonselius’s method, respectively 


(G y ,) exp , apparent shear modulus; (G ? ,) est , corrected shear modulus; 
CV, coefficient of variation 


the correction factor a rapidly reaches a constant value a c 
after a few iterations. The effects of Young’s modulus on a c 
were also investigated numerically. The results (Table 2) 
showed that a c does not vary much for ±25% (±33% for 
LT specimens) changes in Young’s modulus used for FEM 
calculations, especially for RL or TL specimens. 

The shear moduli corrected using a c are shown in Table 
3 together with the apparent shear moduli. It is clear that 
the corrected shear moduli are tolerably equal to each other 
irrespective of the test conditions, whereas this is not the 
case with the apparent shear moduli. Furthermore, it should 
be noted that (G >z ) exp obtained for RL or TL specimens with 
c at 10 mm is closest to (G v ,) est . These results were compared 
with those obtained from the four-point bending tests em¬ 
ploying the same method as we proposed here. They are 
close to each other, as shown in Table 4. Moreover, the 
shear modulus was determined by Fonselius's method 2 in 
which the correction factor is written as 


the above findings, it can be concluded that the newly pro¬ 
posed quasi-simple shear test is useful for estimating the 
shear moduli of wood. Thus, we believe that the correction 
factors approximated from Table 2, which are about 1.1 for 
RL or TL specimens and 1.5 for LR or LT specimens, are 
also used in the experiments conducted on other species. 


Conclusions 

A quasi-simple shear test was undertaken to ascertain 
the feasibility of measuring the shear moduli of wood. The 
shear moduli can be obtained appropriately, provided the 
correction factors a c are determined from a FEM simula¬ 
tion. Approximate values of a c were found to be about 1.1 
and 1.5 for RL or TL specimens and LR or LT specimens, 
respectively. Further work is required to examine the 
validity of this measuring method using other species. 
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